AD-P003  796 


Tutorial  Survey  of  Algoritnas  For  Locating  ana  identifying 
Spatially  Distributed  Sources  and  Receivers 

M.  Morf,  b.  friedlandar  and  J.  Newkirk 
information  Systems  Laboratory 
Stanford  University 
Stanford,  CA  94305 


-T>a  preHSirc  a short  tutorial  survey  of  algo- 
ritnas  for  locating  and  identifying  spatially  dis- 
tributed sources  and  receivers.  The  emphasis  is 
on  aetnods  boat  are  eitner  considered  to  be  very 
basic  or  lend  themselves  potentially  to  distri- 
buted computations,  tne  main  oojective  of  this 
uork.»«e  snail  also  very  briefly  outline  our  own 
approach  to  this  set  of  prooleas. 

X.  Introduction 

in  many  practical  problems  it  is  necessary  to 
determine  the  location  of  signal  (noise)  sources 
frcm  measurements  provided  oy  one  or  more  sensors. 
Typical  applications  include: 

- Acoustic  surveillance  systems  (e.g.,  sonar 
detection  of  iow  flying  aircraft), 

- Seismic  arrays  for  seismic  exploration,  moni- 
toring earthquakes  and  nuclear  explosions,  or 
aetaoting  vehicle  movements, 

- Antenna  arrays  for  radio  astronomy  or  elec- 
tronic surveillance  (e.g.,  direction  find- 
ing), 

- Multiple  raaar  systems  for  detection  and 
tracking. 

The  diversity  of  applications  involving  tne 
target  location  proolem  maxes  a general  unified 
treatment  of  this  subject  quite  difficult.  To 
provice  some  focus  for  our  discussion  we  will  use 
tne  following  sample  problem: 

Consider  a 3mall  numoer  of  sensor  sites 
(perhaps  ten)  dl3triouted  over  a speoiflec 
area.  A number  of  targets  are  present  in  tne 
area  and  tneir  location  is  to  be  estimated 
cased  on  the  data  collected  by  the  sensors. 
The  sensors  measure  signals  wnieh  are  eitner 
aaitted  by  tne  target  (toe  passive  case)  or 
reflected  oy  it  (tne  active  case  which 
requires  target  illumination) . By  processing 
tne  signals  provided  by  the  3ensor,  informa- 
tion about  target  bearing  ana/or  range  can  be 
determined. 

Sometimes  a single  sensor  is  not  capable  of 
measuring  eitner  range  or  bearing,  as  for  example 
witn  omnidirectional  passive  sensors.  Combining 
oata  from  a group  or  array  of  sensors,  however, 
aaxe3  it  possible  to  fine  the  desired  information. 


A sensor  array  of  this  type  may  be  located  at  a 
single  site,  in  which  case  we  conslaer  it  as  one 
unit,  or  it  may  be  distributed  among  many  sites. 

In  otner  cases  the  sensor  sites  can  provide' 
different  types  of  target  related  data,  in  partic- 
ular: 

(1)  Range  only  (ranging  radar,  active  sonar) 

(li)  Bearing  only  (optical,  infrared  sensor, 
direction  finder) 

(iii)  Bearing  and  range  (searon,  tracking  radar) 

(iv)  Target  velocity  (uoppler  radars,  KTI) 

Different  data  types  lead  to  different  location 
estimation  teconiques.  For  example,  range  only  or 
bearing  only  measurements  are  related  to  target 
association  teenniques  (section  2.2).  Bearing  and 
range  data  is  usually  associated  with  tracking 
algor itnms  for  moving  targets. 

Data  from  a single  omnidirectional  passive  sensor 
is  treated  oy  time-of-arrival  mecnods  (seotion 
2.1)  or  beamforming  and  array  processing  teen- 
niques (section  2.4).  The  estimation  metnod  also 
depends  on  the  type  of  signals  provided  by  the 
sensor  site:  conerent/nonconerent,  "raw"  or  fil- 
tered aata  (linear  processing),  data  after  detec- 
tion (nonlinear  processing) , etc. 

Classical  aetbods  of  processing  sensor  data 
have  generally  been  of  tne  centralised  type,  that 
Is,  all  of  tne  sensor  data  was  collected  at  one 
site  and  tnen  processed. . An  alternative  is  to 
process  muon  of  tne  data  at  the  collection  site 
and  to  send  only  tne  relevant  data  to  eitner  a 
central  site  or  (more  generally)  to  the  appropri- 
ate user.  In  Section  3 we  snail  discuss  the  dif- 
ferent types  of  distributed  processing  and  their 
advantages. 

in  the  last  section  we  will  describe  our  own 
approacn  to  tne  csvelopment  of  cistriouted  algo- 
rithms for  the  estimation  of  position,  location 
and  otner  onaracteristios  of  sensors  and  sources, 
(rtorf  et  ai.).  We  will  give  a snort  aesoription 
of  sample  algorithms  of  a fully  distributed  nature 
tnat  nave  oesirabis  features.  »e  snail  also  out- 
line several  of  cue  nonolassicai  approacaes  to 
solving  these  prooleas. 

2.0  Teour.iques  for  estimating  target  location 


fa is  work  wa3  supported  oy  ARPA  under  contract  MLA  This  section  provides  a orief  summary  of  tne 
903-74-0-0179.  solution  teenniques  assoolatea  witn  tne  target 
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location  problem.  Only  the  oasic  ideas  are 
presented;  tne  details  can  oe  found  in  the  refer- 
ences. 


fiae-of- Arrival  sscimation 


A signal  emanating  from  a remote  source  and 
measured  m the  presence  of  noise  at  two  spatially 
separated  sensors  can  be  modeled  as 


x,(t)=s(t;+n,(t) 


x2it)3A3(t+0)+n2(t) , 


wnere  s(t),  n,(c),  n,(t)  are  assumed  to  be  sta- 
tionary, independent,  random  processes.  Cne  com- 
mon aatnou  of  estimating  tne  time  delay,  0,  is  to 
compute  the  cross  correlation  function 


Rx.Xjaiix, (t)Xj(t-W} 


It  follows  directly  tnat 


Rx,x2(r)sAR9S(T>-0) 


wnere  H ( ■ ) is  the  signal  autocorrelation  func- 


tion. An  important  property  of  autocorrelation 


functions  is  tnat  S,«(T)iR  (0).  Thus,  the  peaic 
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of  "..(T-D)  will  occur  at  ttO.  This  provides 


with  a way  for  finding  tne  delay  by  calculating 
tne  estimated  cross-correlation  function.  If  s(t) 
is  a white  noise  source,  R (•tp-QUSOr-D) , and  the 
peax  will  be  snarpiy  defined.  In  general,  R3a( * ) 
will  oe  "spread  out"  whicn  tends  to  broaden  the 
peat;,  matting  it  more  difficult  to  pinpoint  the 
actual  delay,  furthermore,  wnen  multiple  targets 
(ana  multiple  delays)  are  present,  tne  "tails"  of 
tne  autocorrelation  functions  for  aifferent  tar- 
gets will  be  over lay ed  and  more  difficult  to 
separate.  Thus,  it  is  desirable  to  preprocess  tne 
sensor  measurements  x, , x->  so  that  after 
orossoorreiation  3harper  peaxs  will  result,  as  in 
fig.  In  tne  absence  of  measurement  noise  this 
can  oe  done  ay  passing  x,(t),  x~(t)  tnrougn  a 
"wnitening"  filter  for  s(t),  canoe  the  correlation 
of  s is  removed.  wnen  noise  is  present,  tne 
filter  has  to  taxa  into  account  ootn  signal  and 
noise  spectra. 


v ' ; detector 


Ail  of  tne  information  about  the  target  loca- 
tion are  encoded  in  the  relative  tiae-deiays  of 
tne  various  3ensors.  To  3ee  this  consider  tne 


following: 
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Target-sensor  geometry 


wnere  Dj,  is  th*  relative  delay  between  sensor  i 
and  sensor  J,  and  c Is  the  propagation  velocity. 


r?=(xT-xi)2+(y7-yi)2,i  a 1,2,...N. 


it  can  be  shown  iScnmidt]  that  for  h2j,  the  set  of 
equations  (4), (5)  can  be  rewritten  as  a linear  set 
of  equations  for  x™  yTl  wnere  tne  coefficients  are 
known  quantities*’  (i.e.  written  in  terms  of 
°ij  xi  Tbis  3et  of  equations  can  now  be 
soiled' to  determine  the  target  location. 


The  discussion  above  indicates  tnat  cne  way 
of  solving  tne  target  location  problem  is  to  first 
estimate  tne  tiae-of-arrival  delays  and  tnen  to 
compute  the  location  based  on  tne  geometry  of  tne 
problem  [Hahn].  it  is  possible,  of  course,  to 
combine  these  two  steps  and  develop  an  estimator 
aireotiy  for  the  target  coordinates  (x.  y.)  or,  as 
► gis  more  commonly  dons,  for  its  oearing’ana  range. 
This  leads  to  alternative  estimator  structures, 
typically  using  tne  maximum  lixelinood  approach 
[Bangs  and  Scnultneissj , [MacDonald). 


multiple 

correlator 


Target  Association  Techniques 


fig.  A time-delay  estimator. 


Different  cnoices  for  tne  pre-fiiter  are  pos- 
sible depending  upon  toe  performance  criterion 
cnosen  oy  tne  designer:  the  lixelinooo  function 
(hann  and  Tretter],  tne  deflection  function  [Knapp 
ar.d  Carter],  etc.  it  should  be  noted  tnat  several 
estimator  structures  oesides  tr.e  multiplier- 
correlator  estimator  nave  been  saveloped. 


A special  type  of  problem  arises  wnen  multi- 
ple sensors  whicn  measure  range  but  not  asimutn 
(or  vice  versa)  are  used  to  estimate  target  loca- 
tion. if  only  a single  target  is  present,  its 
location  is  found  oy  muitilateration.  for  exam- 
ple, if  azicutn  measurements  from  several  sensors 
are  available,  one  nas  only  to  compute  tne  inter- 
section of  tne  various  lines-of-signt. 


Fig  3:  The  target  association  problem 


£acn  sensor  is  assumed  to  nave  detected  the 
two  targets  I,  However,  tnese  detections  are 
not  properly  associated,  it  is  not  known  wnicn 
measurement  of  each  sensor  corresponds  to  wnicn 
target,  Thu s,  it  is  necessary  to  associate  tar- 
gets uita  sensor  measurements  Before  estimating 
tne  target  locations  (in  fact  tne  association  and 
location  prooleos  are  addressed  simultaneously) . 
:ioce  also  tnat  if  taere  are  more  targets  tnan  sen- 
sors, amoiguities  ("gnost"  targets)  may  result. 
In  fig.  a,  if  was  not  taere,  tne  measurements 
of  a,  and  i2  would  oe  consistent  witn  tne  assump- 
tion taat  tne  targets  are  at  g,  g,  rataer  tnan  at 
1-  - • c 

several  3cnenes  nave  Deen  proposed  to  solve 
tne  target  association  problem,  and  tney  are 
briefly  cescrioed  oeiow. 

',i3t  Forming 

Pick  a pair  of  sensors  and  compute  all  tne  inter- 
sections of  tneir  iine3-of-signt  to  potential 
targets  (i.e.  directions  in  wnicn  tney  detected 
scoetning) . Tnese  intersection  points  are  poten- 
tial target  locations,  sow  picic  a tnira  sensor 
ano  cnecK  wnetner  its  iines-of-signt  pass  througn 
any  of  tne  intersection  points.  If  not,  delete 
tr.e38  points  from  tne  list  of  potential  targets, 
oy  proceeding  tms  way  witn  tne  otner  sensors,  tne 
list  will  finally  include  only  tnose  target  loca- 
tions wnicn  are  consistent  witn  ail  tne  observa- 
tions. it  scbuid  oe  empnasizec  tnat  tnis  is  a 
nigniy  simplified  description  of  more  realistic 
list  forming  aigoritnas. 


Back-projection  or  5pac«-oearea 


Ibe  space  to  be  searcned  is  divided  into  ceils  of 
a size  corresponding  to  tne  system  resolution, 
lae  number  i is  adoed  to  those  ceils  of  tne  space 
wnicn  lie  along  tne  line  of  signt  of  a given  sen- 
sor detection.  This  process  is  repeated  for  ail 
iines-of-signt  of  ail  sensors.  As  can  oe  seen 
from  Fig.  U tne  target  locations  can  oe  identified 
as  tnose  having  tne  nignest  nuaoer  (=  tne  nuaoer 
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of  sensors)  written  m them,  note  tnat  relatively 
higo  values  can  oe  ootained  at  locations  otner 
tnan  tnose  of  tne  real  target  ("gno3t"  -ocations). 


\ , 


Fig.  a:  Association  oy  oaak-projectlon 


Image  reconstruction  techniques 

It  is  possible  to  view  tne  sensor  measurements  as 
line  integrals  tnrough  an  image  consisting  of 
"points  of  light"  at  tne  target  locations.  The 
line  -integrals  are  over  the  infrared  eaisslvity 
map  (for  IB  sensors)  or  tne  radar  reflectivity  map 
(for  ranging  radars).  The  problem  of  reconstruct- 
ing images  from  tneir  line-integral  projections 
nas  been  extensively  treated  in  literature  Le.g., 
Stocks  ana  Qi  Chlro,  Horn),  hecently  it  was  snown 
how  these  techniques  can  oe  applied  to  the  target 
association  problem  (Friedlander  at  al. , Denton  it 
ax.)  oy  reconstructing  the  "brigntness  map"  of  tr.e 
area  under  surveillance  and  identifying  targets  as 
tne  "bright  spots".  It  snouid  be  noted  tnat  the 
image  reconstruction  aetnod  requires  that  tne  sen- 
sors provide  the  actual  energy  measured  in  eaca 
direction  (range)  and  not  Just  target/no  target 
information.  The  detection  takes  place  after  pro- 
cessing tne  information  from  all  tne  sensors;  in 
list  forming  ano  space-searcn,  only  tne  results  of 
the  detection  performed  at  eacn  individual  sensor 
are  passed  on. 

1.3  Spectral  estimation 

Hultiple-sensor  measurements  can  oe  con- 
sidered as  samples  of  a time-space  function  y(t ,£) 
wnere  ? represents  a point  in  3-D  3pace.  The 
notions  of  (temporal)  correlation  function  ar.a 
(temporal)  spectral  density  can  be  extended  to 
time-space  functions  of  tnis  type,  ne  define  a 
random  field  v(t,£)  as  stationary  ana  Homogeneous 
if  styU.it)}  =3  ana 


TAJ 


Mt,t*  ,£,£*}=6{y(t,£Jy(f  ,£•)) 
s"v-t t,z) 


(b) 

wnere 


Thus,  target  tearing  estimation  reduces  to  tne 
profiles  of  estimating  ?(u,ju  from  the  aeasurenents 
where  r.  represents  sensor  locations  ana 
t,  are^the  sampling  times  of  tne  output  of  tnat 
sensor. 


■fct-t',  £*£.-£’ 

Any  Homogeneous  ransom  field  nas  a spectral 
representation 


Many  spectral  estimation  techniques  have  been 
used  in  chis  context;  several  are  described  in  tne 
references. 

2. 4 be am forming  and  Array  Processing 


wnere  £*(*, £=(x,y,z)  and  z(£n£H)  is  a 
random  function  with  certain  properties.  The 
correlation  function  can  oe  represented  by 

y audit  (d) 

— o* 

where  ? is  the  spatial-temporal  spectral  density. 
as  in  tne  temporal  case,  an  inversion  formula 

nolds: 

a(tr ,£}  afrdc.  (9) 

As  an  Important  example,  consider  a mono- 
chromatic (single  frequency)  piano  wave  at  tem- 
poral frequency  ul  propagating  in  tne  direction 
given  by  a unit  vector  ft  at  velocity  a.  Such  a 
wave  is  represented  by 


Pernaps  tne  most  common  operation  in  processing 
signals  in  a senior  array  is  that  of  beamformlng. 
Antenna  arrays  (for  radar,  communication,  etc.) 
and  acoustic  sensor  arrays  (sonar)  are  typical 
examples  of  beamroreing.  Beanforming  consists  of 
a summation  of  time-delayed  (or  pnase-snifteo) 
versions  of  tne  sensor  outputs,  l.e. 

N 

Z(t)a  $ y(t-t  &),  (12) 

1*1  1 ^ 

wnere  t^  represents  the  time-ceiay  for  sensor  i 
ana  it3  location.  Proper  onoice  of  tne  delays 
t,  enhances  the  signals  received  from  a particular 
direction  and  .attenuates  signals  from  otner  direc- 
tions. This  operation  is  tne  spatial  equivalent 
of  a temporal  narrow  handpass  filter. 

The  output  Z(t)  of  the  ceaaforaer  is  a 
(scalar)  time  function,  unlch  is  processed  so  as 
to  obtain  a measure  of  tne  signal  energy  in  an 
optimal  way.  The  moat  coman  processing  scneses 
include 

Matched  filtering 

Wiener  filtering  or  least-squares  estimation 

Maximum  livelihood  estimation 


/<  t.£)*exp(  jutjt+ig'IU 


(10a) 


wnere 


COb) 


it  is  easy  to  verify  that  for  tnis  space-time 
function  we  nave 


.i(^,£)=e 


-jCuijtvjtg'z) 


(’’a) 


and 


P(m,i;)s^(l*-m3,jt-)Jo)  (lip) 

wnicn  is  a delta  function  located  at  temporal  fre- 
quency m3  ana  spatial  frequency  (or  wave  number) 
■Cj.  This  example  indicated  now  ?(u,vj  provides 
information  regarding  tne  direction  (and  velocity) 
of  propagation  of  waves. 


This  linear  filtering  is  often  followed  oy  a 
nonlinear  operation,  e.g.  squaring  and  integra- 
tion. Tne  order  of  linear  filtering  ana 
transforming  (also  a linear  operation)  is  often 
reversed,  tne  exact  structure  depending  on  tne 
application.  The  implementation  of  toese 
processes  is  usually  done  in  tne  rourier  domain 
(with  phase  shifts  replacing  time  delays)  but  time 
domain  implementations  are  also  used,  a sample  of 
the  vast  literature  on  beamforming  ana  tne  associ- 
ated signal  processing  (refened  to  as  "array  pro- 
cessing"! is  given  in  the  bibliograpny. 

A ciass  of  array  processor  of  partisualr 
interest  are  the  different  types  of  adaptive 
arrays.  The  need  for  adaptive  arrays  arises  for 
many  reasons.  Some  examples: 


Hull  steering,  to  minimize  interference  from 
sources  other  taan  tne  target  of  inter»st. 


The  point  of  tnis  discussion  is  tnat  (spa- 
tial) spectral  estimation  is  a way  of  estimating 
target  bearing,  since  if  we  plot  ?(u,£)  in  the  * 
plane  (for  a fixed  «•),  it  will  tend  to  oe  concen- 
trated around  the  point  jt  wnicn  corresponds  to  tne 
direction  of  tne  wave  propagation  and  ner.ce  the 
searing  of  tne  target  (tne  source  of  tnese  waves). 


Adaptive  filtering,  to  mantle  umcnown  noise 
and  signal  statistics. 

Adaptive  beamforming.  Seamforaing  requires 
precise  icnouledge  of  tne  sen-op  locations 
(within  fractions  of  a wavelengta)  m order 


mter- 


tnat  "ae  steering  aelays  be  computed.  Rela- 
tively saall  errors  can  lead  to  serious  per- 
formance degradation.  Thus,  unen  sensor 
locations  are  imprecisely  Known  or  are  con- 
stantly changing,  a fixed  processing  scneae 
is  infeasible. 

;.  The  Heed  for  Distributed  Computation 

There  are  many  advantages  in  distributing 
coaputations  for  a large  sensor  network;  some  of 
tae  main  arguments  are  tne  following: 

( 1 ) Reduction  of  Computational  Complexity 

Distributed  processing  is  often  used  as  a 
means  for  solving  problems  related  to  large-scale 
systems.  This  approch  leads  to  the  decomposition 
of  a nign-dimensional  problem  into  a sequence  of 
snaiier-diaensional  ones.  Tnis  often  results  in 
considerable  computational  savings;  many  fast 
aigoritoos,  sucn  as  tne  Fast  Fourier  Transform, 
are  of  tnis  type.  Also,  certain  iarge-scale  prob- 
lems siapiy  cannot  oe  solved  in  a direct  manner 
(e.g.,  inversion  of  very  large  matrices)  and  ways 
nave  to  ae  found  to  decompose  tne  problem  into 
smaller  parts  tnat  can  oe  nandled.  This  can,  of 
course,  oe  cone  in  a centralized  manner,  but  the 
3i3triouted  approach  often  leads  to  natural  decom- 
positions and  vaiuaole  insights. 

U)  Reliability 

Distributed  systems  have  good  properties 
from  a reliability  standpoint  due  to  their 
icaerent  parallelism.  Faiiuru  of  a computational 
module  does  not  necessarily  result  in  system 
failure  since  toe  computational  load  can  be  re- 
distributed among  the  remaining  modules.  Thus,  a 
distributed  system  may  have  tne  ability  to  recon- 
figure and  continue  operation.  Depending  on  the 
type  of  tne  system  and  its  structure,  its  opera- 
tion after  reconfiguration  may  be  at  a reduced 
performance  level.  (This  would  oe  the  case  If  the 
remaining  computational  resources  were  insuffi- 
cient to  complete  tne  solution  of  tne  problem,  or 
if  tne  loss  of  a computational  module  was  associ- 
ated witn  tne  loss  of  a sensor  site.)  Tne  system 
displays  graceful  degradation  of  performance, 
waicn  contrasts  witn  tne  "catastrophic"  failure 
mode  of  centralized  systems. 

(o)  Flexibility 

The  distributed  nature  of  computations  is 
often  associated  with  distributee  system  struc- 
tures; sucn  a system  structure  aignt  be  a collec- 
tion of  sensor/cooputer/coaounication  modules 
interconnected  In  a network.  This  organization 
leads  to  a ver 7 flexible  structure , possessing 
aesiraeie  properties  wnicn  are  not  always  present 
in  a centralized  system: 


Easy  system  growta 

Tne  capcity  to  candle  topological  onanges  in 
tne  system  structure  (e.g.,  adding  or  delet- 


ing nodes  during  aaintainance  witnout 
rupting  system  operation) . 

The  possibility  of  incorporating  many  combi- 
nations of  resources,  witn  variable  perfor- 
mance levels,  as  determined  oy  tne  needs  of 
each  user. 

- The  combination,  in  a single  network,  of  many 
types  of  information  sensors. 

4.  Our  Approaones 

Distributed  processing  has  by  now  become  a 
term  teat  is  applied  to  many  types  of  systems  and 
is  not  very  weil  defined,  since  tne  distriouted 
sensor  net  prooiem  can  oe  veil  describee,  see  e.g. 
[1SL  DoN  Report;,  we  snail  use  it  as  a basis  for 
defining  distributed  processing.  *e  consider 
three  computational  organizations  which  coula  ce 
used  in  tnis  context:  centralized,  inuepenoently 
distriouted,  and  cooperatively  distributed. 

1 ) Centralized  — all  sensor  data  is  passed 
to  a central  site,  wnere  computation  is  per- 
formed, and  tne  pertinent  resuits  are  taen 
returned  to  tne  appropriate  remote  sites. 

2)  Independent  — ail  sensor  information  is 
communicated  to  every  other  3ite,  and  each 
site  then  makes  its  oest  estimate  of  tne 
environment. 

3)  Cooperative  — sites  exchange  processed 
information,  and  at  most  partial  sensor  data. 

It  is  tae  second  and  third  organizations  tnat  are 
normally  referred  to  as  distributed  organizations, 
and  tne  third,  in  particular,  tnat  we  consider  to 
be  of  greatest  interest. 

The  Eearcn  for  Distributed  Algorithms 


Centralized  algorithms  are  now  quite  veil 
understood,  as  a perusal  of  tne  extensive  litera- 
ture indicates,  nowever,  very  few  attempts  nave 
been  made  to  unify  and  integrate  ail  tuese  dif- 
ferent approaches  ana  results.  A typical  cook  on 
radar  or  sonar  signal  processing  is  a ratner  aa 
hoc  collection  of  data,  oecnods  ana  cneory  (remin- 
iscent of  a cook-oook) . To  an  outsider  of  this 
field  it  is  extremely  difficult  to  get  a conerent 
picture  and  to  maxe  intelligent  cnoices  :n  apply- 
ing tnese  methods  in  the  design  of  systems.  A 
systematic  representation  of  this  xnowlege,  oy 
itseif  a tremendous  task,  is  required  ir.  oruer  to 
make  effective  use  of  tne  availaoie  alternatives. 
It  is  very  tempting  to  suggest  the  development  of 
an  "expert  support  system"  combining  and  extending 
recent  approaches  in  Al,  cata-oase  management  ir.c 
related  fields. 

Our  approach  to  tne  development  of  (coopera- 
tively) distributed  algcritnms  can  se  summarized 
under  tne  following  headings: 

'■  Partitioning  of  optimal  centraiicec 
algorithms,  sucn  as  Kaximum-i.ikelir.coc, 
Extended  Salman  Filtering,  and  seam 
Forming.  This  approach  is  useful  wner. 
tr.e  system  under  consideration  can  co 


divided  Into  subsystems  vita  sparse 
Interaction. 

2.  Application  and  extension  of  methods 
developed  in  tae  context  of  decentral- 
ized control  and  estimation,  e.g.  team 
decision  theory  and  differential  games, 
Hierarchical  and  multilevel  systems, 
aggregation  aetnods,  singular  perturba- 
tion and  otner  perturbation  techniques, 
periodic  coordination  and  spatial 
dynamic  programming  (sdp). 

3.  Development  of  new  optimal  distributed 
algoritms  for  specific  subsets  of  sensor 
data.  For  example, 

Time/Frequency  Difference  Of 
arrival  (TuGA/FDQA)  data,  using 
ARM*  modeling. 

Range  only  cr  Angle  only  data, 
using  image  reconstruction  tech- 
niques or  distributed  versions  of 
the  bacxprojection  technique 
descrioed  in  Section  2.2. 

- Range  and  Angie  data,  using  non- 
linear intimation  techniques. 

Mixed,  possibly 

inconsistent/ incomplete  data,  using 
a hierarchical  approacn. 

4.  Advanced  Concepts 

The  physical  proolem  of  locating 
and  identifying  sources  aas  much 
natneaatical  structure;  for  exam- 
ple, it  is  heavily  dependent  upon 
tae  o noice  of  coordinate  systems. 
nan  euclidean  Geometry  and  Non 
Classical  Statistics  (Non-Gaussian) 
vui  very  prooabiy  oe  of  great 
oeneflt. 

- Signal  processing  of  one- 
dimensional  signals  is  a very  weli 
developed  field;  spatial  and  other 
multi-  dimensional  problems,  how- 
ever, require  more  advanced 
mathematical  tools. 

- Our  preliminary  investigations 
indicate  that  most  candidate  algo- 
rithms for  distributed  processing 
require  high  communication 
sanawidtns.  As  an  alternative 
approach,  we  are  Investigating  Pro- 
babilistic Algorithms;  these  algo- 
rithms potentially  require  lower 
bar.dwidths,  are  naturally  suited  to 
parallel  and  distributed  organiza- 
tions, and  they  can  os  very  robust. 

- from  a systems-desiga  perspective 
one  snouia  consider  interactions 
between  software  and  hardware 
architecture  early  in  the  cevelop- 


nent  of  algorithms.  For  this  rea- 
son, we  are  considering  the  poten- 
tial impact  of  VLSI/VHSl  designs. 

Using  these  approaches,  examples  of  fully 
distributed  processing  and  communication  algo- 
rithms can  be  proposed.  One  suen  example  is  tne 
combination  of  the  TDOA  approacn  [Scnaldt],  the 
distributed  estimation  algorithm  in  [ iSL-DSN 
Report]  and  a distributed  protocol  a la  [Merlin 
and  degall J . These  algorithms  nave  the  desired 
robustness  and  low  communication  banduidtbs  taat 
characterize  desirable  distributed  algorithms. 
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